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Abstract: We formulate and demonstrate experimentally the high- 
resolution spectral method based on Bloch-wave symmetry properties for 
extracting mode dispersion in periodic waveguides from measurements of 
near-field profiles. We characterize both the propagating and evanescent 
modes, and also determine the amplitudes of forward and backward waves 
in different waveguide configurations, with the estimated accuracy of several 
percent or less. Whereas the commonly employed spatial Fourier-transform 
(SFT) analysis provides the wavenumber resolution which is limited by the 
inverse length of the waveguide, we achieve precise dispersion extraction 
even for compact photonic structures. 
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1. Introduction 



The latest advances in near-field measurements open new opportunities for experimental char- 
acterization of pulse dynamics inside complex photonic structures |lj]. The correlation measure- 
ments of the near-field probe signal and the reference pulse enable the determination of both 
amplitude and phase of the electric field at distinct time frames. These data enable direct visu- 
alization of pulse propagation in real space. Additionally, it was shown that spectra of the field 
profiles obtained through spatial Fourier-transform (SFT) can be used to reveal the dispersion 
of modes in optical waveguides |2|. Such analysis is especially valuable for photonic-crystal 
waveguides, where dispersion curves can possess complex features such as turning, inflec- 
tion, and anti-crossing points OS). However, there exists a fundamental limitation on results 
obtained with SFT: Ak > 2n/L, where Ak is the resolution of the wavenumber, and L is the 
structure length. This severely restricts the potential to determine the mode dispersion. 

High-resolution spectral methods have been developed for analysis of temporal dynam- 
ics J5]|6l, overcoming the limitations of Fourier-transform method. An implementation of such 
approach for dispersion extraction of photonic-crystal waveguide modes was recently shown 
with numerical simulations 0. In this work, we reveal that the spatial spectral analysis can be 
enhanced through the application of generic Bloch-waves symmetry properties in arbitrarily- 
shaped periodic dielectric waveguides. In particular, we show that the required structure size 
for dispersion extraction can be further reduced, offering additional advantages for characteri- 
zation of short waveguide sections, which may form parts of compact photonic-crystal circuits. 
Furthermore, we demonstrate accurate extraction of real and complex wavenumber values from 
near-field experimental measurements, precisely determining the mode dispersion both inside 
the transmission band and within the photonic band-gap regions. 

The paper is organized as follows. In Sec. [2] we formulate the general approach for single- 
and multi-mode periodic waveguides based on Bloch-wave symmetries. Then, in Sec. [3] we 
demonstrate dispersion extraction for a single-mode periodic waveguide and transmission line, 
and in Sec. [4] present an example for two-mode homogeneous waveguide. The examples in 
Sees. [3] and |4] are based on experimental measurements at the microwave frequencies, confirm- 
ing the high accuracy and robustness of extraction procedure under the presence of noise in 
experimental data. Conclusions and outlook are presented in Sec. [5] 

2. Dispersion extraction based on the Bloch-wave symmetries 

The goal of high-resolution spectral analysis J5) is to determine the frequencies or wavenum- 
bers of a finite total number of modes (M) which primarily determine the system evolu- 
tion. Let us consider how this methodology can be applied to extract dispersion for a peri- 
odic waveguide which supports M modes in particular frequency range. Since each of the 
modes of a periodic waveguide satisfies the Bloch theorem HQ, the complex electric field 
envelope of a waveguide mode with the index m at the frequency co can be expressed as 
E(x,y,z;co) = y/ m (x,y,z',(o)exp(ik m z/d). Here k m is the Bloch wavenumber, x and y are the 
orthogonal directions transverse to the waveguide, z is the direction of periodicity, d is the 
waveguide period, and i//,„ is the periodic Bloch-wave envelope function: y/„,(z) = y/ m (z + d). 
Then, the total field inside the waveguide can be presented as a linear superposition of propagat- 
ing modes and radiative waves: E(x,y,z; co) = Lm=i a m \j/ m (x,y,z; co) exp(ik m z/d)+w(x,y,z; co), 
where a m are the mode amplitudes and w(x,y,z', co) is the radiation field due to the excitation of 
non-guided waves. 

The key question is how to use the general spectral methods J3] to determine the mode 
amplitudes a m and wavenumbers k m which provide the best fitting for the experimental 
measurements of the electric field distribution. One approach is to use the periodicity prop- 
erty of Bloch-wave envelopes and represent them as infinite Fourier series, \j/ m (x,y,z',d)) = 



LitTloo Wm,s{x,y,C0)exp(i2nsz/d). Then, by performing the measurements along a line with 
fixed transverse positions (xo,yo), me problem is reduced to finding harmonics with the ampli- 
tudes yf ms (xo,yo;Ob) and the corresponding wavenumbers k m>s = (k m + 2ns). In case of Bloch- 
waves with smooth profiles, terms with s ^ are small allowing efficient dispersion extraction. 
It appears that the latter condition was satisfied in the analysis of Ref. Q, enabling accurate 
dispersion calculation. However, if the Bloch-wave profiles vary strongly over the unit cell, 
their spectrum would span a broad range of s, substantially complicating the spectral analysis. 

We suggest that, in the general case of arbitrarily-shaped periodic waveguides with complex 
Bloch-wave profiles, accurate dispersion extraction can be performed using the electric field 
measurements at a set of N periodic locations: U„ = E(xo,yo,Zo + n d; co), n = 1 : N. The posi- 
tions (xo,yo) can be chosen at the locations with the maximum field amplitudes, reducing the 
effect of noise in measurements. Taking into account the periodicity property of Bloch-wave 
envelopes \j/ m , we introduce new variables A m = a m \j/ m (xo,yo,Zo', co)exp(ik m zo/d), and obtain a 
set of equations: 

M 

U„ = £ A m exp(ik,„n) +w n . (1) 

m=l 

Then, considering the radiation waves to be small, we can determine the Bloch-wave parame- 
ters which describe most accurately the measured field profile using established mathematical 
algorithms [5 9|. In general, the number of measurements shall be equal to or exceed the num- 
ber of unknowns, N > 2M. 

We note that the spatial spectral analysis can benefit through the application of additional 
constraint due to the symmetry of forward and backward modes in dielectric waveguides, which 
wavenumbers are related as k%m = — &2m-i El - Then, the required structure size for dispersion 
extraction is reduced: N > 3M/2. In order to take advantage of this relation, we use the least- 
squares method in place of general spectral algorithms [5 9|. Specifically, we seek the values of 
A m and £,„ which minimize the functional W = L^=i \ w n\ 2 / Y,n=i \U n \ 2 - For given wavenumbers, 
the minimum Wa({A:,„}) = min^W is achieved when dW/dA,„ = dW/dA* n = 0. It follows 
that the optimal amplitudes satisfy the linear matrix equation C A = B, where components of 
vector A are the optimal amplitude values, components of the matrix C are C pq = exp[i(k p — 
k*)n], and vector B components are B p = L n ^« ex P( — ik*n) for p,q = 1 : M. We can show 
that WA({k m }) = W A= ^ = 1 — L/;^/>.Bp/Ln \U„\ 2 . The remaining task is to find the absolute 
minimum W m ; n = tninj^ Wa (note that, by definition, Wa is real and positive), and this can be 
done numerically, for example by using the 'fminsearch' function in Matlab. 

Due to the possible presence of radiation and noise in experimental measurements, the opti- 
mal value of the functional W m ; n will be larger than zero. Most importantly, this value can be 
used to estimate the accuracy of the dispersion extraction. Specifically, we define the confidence 
interval of wavenumbers as those corresponding to WA(k m ) < 2W m i n . 

3. Single-mode periodic waveguide 

We now illustrate the application of the general approach for dispersion extraction in the 
regime when a periodic dielectric waveguide is single-moded. Then, we take M = 2, ac- 
counting for one forward and one backward mode, which propagation constants are related 
as k2 = —k\. In this case, the required number of measurements at periodic locations for dis- 
persion extraction is > 3. For three measurements (N = 3) we can find an analytical so- 
lution to W = 0: the wavenumbers are k\ = —ki = cos~'[([/i + U3) / (2U2)], and the corre- 
sponding amplitudes of counter-propagating waves are A; = [2U2 + i{U\ — U3)/ sin (k)] /4 and 
Ao = [2I/2 + 1(1/3 —U\)/ sin(fc)]/4. The larger number of measurements allows one to esti- 
mate the accuracy due to the noise and radiation affecting experimental data. It is instructive 




Fig. 1. (a) Sketch of the periodic waveguide made of dielectric pillars, (b) Dispersion re- 
trieved from experimental measurements, shown are both the real and imaginary parts as 
indicated by labels. Grey shading marks the confidence interval, (c) Retrieved amplitudes 
of experimentally excited forward (A+) and backward (A_) waves, (d) Numerically calcu- 
lated dispersion, (e) Normalized SFT of experimentally measured electric field profile at 
rod centers along the waveguide. 



to compare this procedure with the SFT, where the location of the spectral peaks can only be 
determined with accuracy no better than Ak > 2n/N, and for N = 3 practically no information 
on wavenumber values can be obtained. 

We now demonstrate dispersion extraction from experimental data. The measurements have 
been performed for microwave frequencies, where the amplitude and phase profiles of electric 
field distributions excited by vector network analyzer were sampled by a monopole near-field 
probe connected to the second port of the vector network analyzer. 

First, we study wave propagation through an array of 20 dielectric pillars, as schematically 
shown in Fig.[TJa). This system can be used as a testbed for the key features of wave propaga- 
tion in periodic waveguides due to the presence of photonic band-gaps and associated changes 
in dispersion characteristics iflOl . Each pillar is 1cm high, has a diameter of 4.18mm with the 
refractive index of 2.5, and they are sandwiched between the metal plates in the vertical direc- 
tion. The field measurements are performed at each of the rod centers. The extracted dispersion 
for TM polarization is presented in Fig. QIb). We see that the essential features of the pho- 
tonic band-gap are fully recovered. As the frequency is increased, the wavenumber approaches 
the edge of the Brillouin zone. Then, the imaginary part increases, indicating the exponential 
decay of the wave due to the photonic band-gap. Note that, close to the gap-edge, the real 
and imaginary parts of the wavenumber demonstrate the square-root dependence on frequency 



Fig. 2. (a) Photograph of transmission line, (b) Dispersion retrieved from experimental 
measurements (dashed lines) and calculated numerically (solid lines), (c) Retrieved ampli- 
tudes of experimentally excited modes with positive (A+) and negative (A_) wavenumbers. 
(d) SFT of experimentally measured electric field profile along the line. 

detuning, in agreement with the general properties of periodic photonic structures [8|. Addi- 
tionally, we also extract the amplitudes of Bloch waves which are excited experimentally, see 
Fig. 0Jc)- The nontrivial dependence of amplitudes on the frequency is due to the reflections 
at the structure boundaries. It is interesting to see that only one Bloch wave is present for a 
range of frequencies inside the band-gap, since the mode decays quickly and cannot excite the 
other evanescent wave at the opposite end of the waveguide. We also calculate numerically 
the modal dispersion considering the two-dimensional photonic structure with the effective re- 
fractive index of 2.4296, by computing eigenmodes of Maxwell's equations using the freely 
available MIT Photonic Bands (MPB) software package utilizing the plane-wave expansion 
method ifTTl . see Fig.[TJd). Our experimentally extracted dispersion curves perfectly agree with 
the calculated dependencies. For comparison, we also present the SFT spectrum of the electric 
field measurements in Fig. QI e )- This spectrum can only indicate the general dispersion trend 
with very rough wavenumber resolution, and only inside the transmission band. 

Now we analyze the dispersion properties of the left-handed transmission line shown in 
Fig. |2 a), which characteristics and nonlinear response were previously considered in Ref. lfT2l . 
It was shown that for this structure, co(k) = co(—k) similar to dielectric waveguides, and there- 
fore we use this symmetry in dispersion extraction. The line has 20 periods to ensure that 
end effects do not dominate the response, and the field distribution is registered at periodic 
locations along the line. We observe an excellent performance of extraction method for this 
periodic structure, as demonstrated by close agreement of dispersion curves based on experi- 
mental data and calculated numerically, see Fig. |2jb). We also determine the mode amplitudes 
and find that the dominant propagating mode has negative wavenumber at lower frequencies, 
and positive wavenumber at higher frequencies above the stop-band. This result agrees with the 
nature of modes in the transmission line lfl2l . As in the previous example, we find that very 
limited information on mode dispersion can be obtained with SFT, shown in Fig|2jc). 

4. Multi-mode waveguide 

We now demonstrate the dispersion extraction for a multi-mode waveguide. Specifically, we 
consider a homogeneous (non-periodic) perspex dielectric slab waveguide, shown schemati- 
cally in Fig. |3 a). The waveguide is 3cm wide and 20cm long, and it is sandwiched between the 




Fig. 3. (a) Sketch of perspex dielectric slab waveguide, (b) Dispersion retrieved from ex- 
perimental measurements (dashed lines) and calculated numerically (solid lines). Shaded 
are the regions above the free-space light line, and below dielectric light line, where guided 
modes do not exist, (c) SFT of experimentally measured electric field profile along the 
waveguide. 



metal plates in the vertical direction. Under such conditions, the dispersion of TE modes can be 
calculated using expressions for planar waveguides lfl3l with the effective dielectric constant, 
which value for our system is found to be £eff = 2.1. In Fig. [3jb) we show the numerically 
calculated dispersion for the first two modes with solid lines, presented for both forward- and 
backward-propagating waves with positive and negative wavenumbers, respectively. 

In experiment, we measure the electric field amplitude at N = 71 points with d = 2mm spac- 
ing along the waveguide. The transverse antenna position is shifted away from the waveguide 
center by xq = 1cm, as indicated by dashed line in Fig.(3ja). This position is chosen since at the 
waveguide center the amplitude of the second mode vanishes as its profile is anti-symmetric. 
Then, we process the data by running the retrieval algorithm for M = 4 total number of modes, 
and we analyze the domain of real-valued wavenumbers kj in order to characterize the guided 
waves. The results of dispersion extraction are shown with dashed lines in Fig.[3jb), revealing 
excellent agreement with numerical calculations. Indeed, we observe that for all frequencies 
the relative error W m \ D remains below 8%. 

For comparison, in Fig. 0c) we show the SFT of electric field amplitudes. The spectral 
resolution of this plot is severely limited due to the small waveguide length. Additionally, the 
amplitude of the second excited mode is much smaller than the fundamental mode, and even 
the existence of the second mode cannot be reliably determined from the Fourier spectra 

5. Conclusion and outlook 

In conclusion, we have presented a generic approach based on high-resolution spectral methods 
combined with specific Bloch-wave symmetries for the extraction of mode dispersion and am- 
plitudes using near-field measurements in periodic waveguides. We have demonstrated its suc- 
cessful application through processing of experimental measurements for a variety of waveg- 
uide configurations, observing very good accuracy and high robustness to experimental noise. 
One interesting problem for the future studies is to explore the possibility to extend this ap- 
proach for two-dimensional Bloch waves. 



